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Abstract The Odorrana schmackeri sensu lato is a group 
of cascade frogs that are widely distributed in southern 
China. However, their taxonomy has long been unresolved 
and several cryptic species have been reported in the past 
decade. Previous analyses revealed the existence of three 
sympatric species (1.е., О. schmackeri, O. hejiangensis and 
O. sp2) at the type locality of O. schmackeri (Gaojiayan 
Town, Yichang City, Hubei Province, China). However, 
the lack of an original description for the holotype 
specimen (SMF6241) of O. schmackeri poses a confusing 
and controversial problem in species nomination within 
this group. Here, we reevaluate the taxonomic affinities 
of SMF6241 and examine the distinctness of О. sp2 based 
on specimens from the type locality. We use integrative 
approaches that combine morphology with 10 microsatellite 
loci and bioacoustics. Canonical discriminant analysis 
of morphological characteristics revealed three distinct 
groupings (ie, О. schmackeri, O. hejiangensis and О. sp2) and 
attributed SMF6241 to O. schmackeri. Genetic clustering 
analysis of nuclear DNA also identified three corresponding 
clusters. In addition, analysis of interspecific variation in 
morphometrics and acoustics showed divergence between 
O. sp2 and its sympatric congeners. All lines of evidence 
from morphological, molecular and bioacoustic approaches 
suggest the species status of O. sp2, namely Odorrana 
ichangensis Chen, sp. nov. Overall, this study resolves the 
taxonomic assignment of the holotype of O. schmackeri 
and provides an example for species delimitation using 
integrative approaches. 
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1. Introduction 


The two fundamental tasks of taxonomy are species 
delimitation and classification (Sauer and Hausdorf, 2011; 
Boury-Esnault et al., 2013). The odorous frogs belong to the 
genus Odorrana and they inhabit montane streams throughout 
subtropical or tropical regions of East and Southeast Asia, 
ranging across southern China, Japan, northeastern India and 
Southeast Asian countries (Fei et al., 2012; Frost, 2020). This 
genus represents a highly diverse species group containing 58 
recognized species, and more than 20 new species have been 
described in the past decade (Pham et al., 2016; Frost, 2020). 
However, the phylogenetic relationships and species validity 
within this genus have long been controversial due to limited 
and incomplete original description, insufficient type material 
and lost holotype specimens (Hallowell, 1861; Boettger, 1892; 
Matsui, 1994; Bain et al., 2003; Kurabayashi et al., 2010; Chen et 
al, 2013; Rujirawan et al, 2018; Chen et al, 2020). 

Odorrana schmackeri (Boettger, 1892) was erected based on only 
one male specimen SMF6241 (= 1054, 2a) from Gaojiayan Town 
(= Kao-cha-hien) of Yichang (Ichang) City, Hubei Province in 
central China (Boettger, 1892). However, the lack of original 
descriptions, diagnostic characters, pictures and etymological 
meaning of the type specimen limited our knowledge of this 
species for many years (Supplementary Figure S1). Liu and Hu 
(1961) provided a relatively comprehensive description of O. 
schmackeri based on 9 specimens from the area near the type 
locality (ie, Wushan and Xiushan County in Chongqing City), 
and it was accepted universally. Due to similar morphology 
and color patterns, the O. schmackeri sensu lato was formerly 
regarded as a single, widespread species in southern China over 
a long time (Fei et al, 1990; Fei et al, 2012). However, the results 
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of comparative morphological and molecular studies revealed 
substantial species-level morphological and genetic variation 
among geographic populations, suggesting that O. schmackeri 
was paraphyletic with several discrete species, forming "the 
O. schmackeri species complex" (Chen et al., 2013; Zhu, 2016; 
He et al, 2017; Li et al, 2017). Four other distinctive species, ie, 
O. hejiangensis, O. huanggangensis, O. tianmuii and O. kweichowensis, 
have been described in this complex (Deng and Yu, 1992; Chen 
et al, 2010a; Chen et al, 2010b; Li et al., 2018). 

Using 1021 samples from 85 sites across 15 provinces (or 
municipalities) of China, a mtDNA-based phylogeographic 
study (Zhu, 2016) demonstrated that the O. schmackeri sensu lato 
was paraphyletic and included six independent clades (Clade 
А-В): О. schmackeri (Clade A), О. kweichowensis (Clade B; named 
as “О. spl” by Zhu), О. hejiangensis (Clade C), О. huanggangensis 
(Clade E), O. tianmuii (Clade F) and one cryptic species (Clade 
D; ie, O. sp2). The mtDNA topology suggested the monophyly 
of all samples of O. sp2 and a sister-group relationship to O. 
hejiangensis with strong support. The sister-group relationship 
between O. schmackeri and O. kweichowensis and between O. 
huanggangensis and O. tianmuii were also well supported. The fine 
geographical distribution ranges of the species in this complex 
were also revealed by Zhu (2016) and several other studies (He 
et al, 2017; Li et al, 2017; Wang et al, 2019), and the results were 
different from previous studies (Fei et al, 2009; Fei et al, 2012). 
However, due to the lack of morphological, acoustic and nuclear 
evidence, the validity of the cryptic species O. sp2 remains 
unresolved. In addition, field surveys indicated that three species 
of this complex (ie, О. schmackeri, O. hejiangensis and О. sp2) 
were sympatrically distributed at the type locality (Gaojiayan 
Town, Hubei Province) of О. schmackeri (Figure 1). Furthermore, 
the affinity of the holotype (SMF6241) of O. schmackeri to which 
of the three sympatric species is still unknown. 

Here, we used morphological data and performed canonical 
discriminant analysis to address the taxonomic attribution 
of the holotype (SMF6241) of O. schmackeri. We also used 
integrative data of morphological characteristics, ten nuclear 
loci and male advertisement calls to test the taxonomic validity 
of the cryptic species O. sp2. Furthermore, we provided the 
taxonomic descriptions of the new species. 


2. Materials and methods 


24. Sampling Specimens and tissues used in this study were 
collected from field surveys, and the majority of them were kept 
in the Amphibian and Reptile Laboratory of Henan Normal 
University. A total of 190 adult specimens were collected from 
eight localities (2-144 samples per site; Supplementary Table 
51). Specifically, specimens of О. sp2 (51 6 , 31 ? ) were from the 
eight sites, while specimens of О. schmackeri (31 5, 23 9) and О. 
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Figure 1 Photos of the sympatric О. schmackeri, О. hejiangensis 
and O. sp2 and their habitat at Gaojiayan Town, Changyang 
County (Co.), Hubei Province (Prov.). (A) Male 
(HNNU18071139, the holotype specimen) and female 
(HNNU18071155, the allotype specimen) О. ichangensis sp. 
nov. in amplexus. (B) Male (HNNUO090811386) and female 
(HNNUO?9081I385) O. hejiangensis in amplexus. (C) Male 
(HNNU18071132) and female (HNNU18071131) О. schmackeri 
in amplexus. (D) The habitat at Gaojiayan Town. 


hejiangensis (33 6 , 21 9 ) were all from Gaojiayan Town, Hubei 
Province (ie, the type locality). All the samples were measured 
for morphological analysis, and 92 of them were used for 


molecular analysis. 


2.2. Morphological measurements All measurements were 
made with a digital caliper to the nearest 0.1 mm precision from 
preserved adult specimens. The following 24 morphometric 
measurements were taken: snout-vent length (SVL; length from 
tip of snout to vent); head length (HL; distance between angle of 
jaws and snout-tip); head width (HW; width of head at level of 
maxilla articulation); snout length (SL; length from the front of 
eye to Ше tip of snout); nostril-snout tip distance (NSD; distance 
from anterior margin of nostril to snout tip); nostril-eye length 
(NEL; distance between anterior corner of eye and nostril 
center) ; internarial distance (IND; distance between central 
points of nostrils); interorbital distance (IOD; shortest distance 
between medial edges of eyeballs in dorsal view); IFE (distance 
between anterior corner of eyes); IAE (distance between 
posterior corner of eyes); eye diameter (ED; eye diameter from 
anterior to posterior corner); maximum width of upper eyelid 
(UEW; maximum distance between medial edge of eyeball 
and lateral edge of upper eyelid); tympanum diameter (TD; 
maximum diameter of tympanum); tympanum-eye distance 
(TED; distance from anterior edge of tympanum to posterior 
corner of the eye); length of lower arm and hand (LAHL); 
lower arm diameter (LAD; maximum dia meter of lower arm); 
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hand length (HAL; length from elbow to the tip of third finger); 
length of hind leg (HLL; length from vent to tip of the longest 
toe); tibia length (TL; length from outer edge of flexed knee to 
heel ); tibia width (TW; maximum width of tibia ); foot and 
tibiotarsus length (FTL; length from tibiotarsal joint to end 
of fourth toe); foot length (FL; length from the base of inner 
metatarsal tubercle to the tip of fourth toe); outer metatarsal 
tubercle length (OMTL; maximum length of outer metatarsal 
tubercle); inner metatarsal tubercle width (IMTW; maximum 
width of inner metatarsal tubercle) (see Fei et al., 2009). We 
measured the above 24 morphological variables of collected 
specimens of O. schmackeri, O. hejiangensis and O. sp2. We also 
obtained 18 variables (SVL, HL, HW, SL, NSD, NEL, IOD, ED, 
UEW, TD, TED, LAHL, HAL, TL, TW, FTL, FL and OMTL) 
of the holotype specimen (SMF6241) preserved at the Muséum 
national d' Histoire naturelle, Paris. For the holotype of the 
new species (HNNU18071139), another two tra its (FPW, finger 
disk width at the widest part of the pad of finger Ш and DPW, 
width of distal phalanx of finger III) were also measured in 
addition to the above 24. For comparisons, morphological data 
for O. schmackeri and O. hejiangensis were based on specimens 
in this study and the original descriptions from Liu and Hu 
(1961) and Deng and Yu (1992). The morphological data of O. 
kweichowensis, O. tianmuii and O. huanggangensis were taken 
from Li et al. (2018) and Chen et al. (201 a, b). 


2.3. Microsatellite loci genotyping and genetic structure 
Ninety-two individuals from Gaojiayan Town were used in 
the molecular analysis (48 O. schmackeri, 21 O. hejiangensis and 
23 O. sp2). Total genomic DNA was extracted using a standard 
phenol/chloroform procedure (Sambrook and Russell, 2001). 
A panel of ten microsatellite markers (Table 1) was used (Gao 
et al, 2013; Qiao et al, 2015). Genotypes for the microsatellite 
markers were conducted following Qiao et al. (2015). Micro- 
Checker v.2.2.3 (Van Oosterhout et al., 2004) was used to 
remove genotyping errors as the presence of null alleles for 
each microsatellite locus. Gene pop on web 4.0 was used to 
test Hardy-Weinberg equilibrium and linkage disequilibrium 
between each pair of loci (Raymond and Rousset, 1995). 
Observed (Но) heterozygosity, expected (Hp) heterozygosity, 
the polymorphism information content (PIC), and number of 
alleles per microsatellite locus were calculated using Cervus 2.0 
(Kalinowski et al., 2007). 

We used the Bayesian-based genetic clustering program 
STRUCTURE 2.3 to infer the number of clusters and 
estimate admixture among these sympatric lineages based on 
microsatellite data (Pritchard et al., 2000). A correlated allele 
frequency model allowing for admixture of individuals was 
assumed. Independent runs were performed by setting the 
number of clusters (K) from 1 to 10. Each run comprised 200 000 
burn-in iterations, followed by 1 200 000 generations of the 
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Monte Carlo Markov Chain (MCMC) to estimate LnP (D). The 
successive changes of LnP (D) values and the delta K value were 
used to assess the optimal K value (Evanno et al, 2005; Earl and 
Vonholdt, 2012). 


2.4. Acoustic analysis As a key role in reproductive isolation 
(Hoskin and Higgie, 2010; Wilkins et al., 2013), bioacoustic 
analysis is currently recommended as a powerful tool for 
anuran species delimitation (Lacerda et al., 2011; Alstróm and 
Ranft, 2013; Köhler et al., 2017). We recorded advertisement 
calls of two О. schmackeri and five О. sp2 separately on July 21 
and 26, 2018 at Gaojiayan Town, Hubei Province. Calls of one 
individual of O. hejiangensis were recorded on July 8, 2018 in 
Wanyuan City, Sichuan Province. The local air temperature 
and humidity were measured with a digital thermo-hygrometer 
to the nearest 0.1°C. 

Male calls were recorded with a Marantz PMD661MKIII 
digital recorder connected to a Sennheiser ME 66 
directional microphone capsule and a Sennheiser K6 powering 
module in wave format files with a sampling rate of 44.1 
kHz and 16-bit resolution. The microphone was placed 
approximately 05-1 m away from the calling frogs (Köhler et 
al, 2017). 

Calls were analyzed using the sound analysis software 
Raven Lite 2.0 (Cornell Laboratory of Ornithology, Ithaca, 
NY, USA) and Adobe Audition 3.0 (Adobe Systems, San Jose, 
CA, USA). Spectral and temporal variables were measured in 
a single call. Spectral measurements were taken through fast 
Fourier transform set at a length of 512 points. Sound figures 
were obtained using PRAAT software (Boersma, 2001). 


2.5. Statistical analysis A canonical discriminant analysis 
(CDA) was performed to discriminate species and predict group 
membership of SMF6241, as in other cases (Parsons and Jones, 
2000; Polaszek et al, 2010; Xiong et al., 2015). We first carried 
out the CDA based on 24 morphological variables of measured 
individuals in males and females, respectively. Then, we 
performed the CDA again on the 18 variables from SMF6241 
and all male individuals used above. In the CDA, canonical 
functions can be derived to compute the squared Mahalanobis 
distance between individuals and group centroids (Cerqueira 
and Lemos, 2000). The classification of SMF6241 was then 
allowed to be determined by the smallest squared Mahalanobis 
distance. The statistical significance of discrimination was 
determined by Wilk's test. Figures were plotted using the 
ggplot2 package in R version 3.22 (Wickham, 2009). 

To examine the morphological intraspecific variance and 
screen for useful variables, we ran one-way analysis of variance 
(ANOVA) and pair-wise multiple comparisons by Dunns test. 
Size differences between males and females were tested by the 
analysis of covariance (ANCOVA) with SVL as covariate. Data 
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Table 1 Characterization of 10 microsatellite markers for O. schmackeri, O. hejiangensis and O. sp2. 

Locus Repeat motif Primer sequence (5'-3) (Е, forward; В, reverse) Size range (bp) N, Та/°С Но Hg PIC 
OG-13 (САТА) » Е т 366-430 11 52 0.461 0.752 0.73 
OG-19 (САТА) d nn 233-271 10 54 0.247 0.668 0.668 
ОС-30 (AGAT), | с 272-348 18 56 0.64 0.893 0.893 
OG-35 (AG); о о у is 315-319 3 58 0.18 0.553 0.486 
OG-50 (GRANT Pp CH MEM Rd 196-244 14 56 0.618 0.867 0.848 
OG-62 (АСАТ), : en 154-190 11 56 0.685 0.849 0.826 
Om-14 (ААС), 1 EE = = 147-159 5 60 0.146 0.635 0.583 
Om-54 (AAT), E f 252-261 4 52 0.292 0.64 0.568 
212884 (TAC), | ее 179-191 5 54 0.483 0.548 0.488 
177957 (АААТ), FF 208-220 4 50 0.258 0.464 0.415 


R: ATAAGGGCCCTCTCAGGGTA 


were log-transformed for these statistical analyses. However, 
mean values (+ standard deviation) and ranges of variations 
are presented as non-transformed data. Inter-species differences 
in acoustic parameters were evaluated using a Kruskal-Wallis 
ANOVA on Ranks adjusted by Bonferroni-Holm. All statistical 
analyses were performed using SPSS software (Version 25.0, 
SPSS Inc, Chicago, IL, USA). 


3. Results 


3.1. Morphology The canonical discriminant function 
based on 24 morphometric characters significantly (all P < 
0.05; Supplementary Table S2) differentiated all individuals 
(with a priori species determination) into a specific group (ie, 
100% accurate classification rate in raw data), corresponding 
to О. schmackeri, О. hejiangensis or О. sp2 for both males and 
females (Figures 2A, B). The leave-one-out cross-validation 
tests were further performed to evaluate the accuracy, and the 
classification accuracy was 91.7% in males and 77.6% in females, 
respectively. The CANI (explaining 87.9% and 77.1% of the 
whole variance in males and females, respectively) suggested the 
differentiation of О. hejiangensis from О. schmackeri and Озр2, 
mainly contributed by TL and TD in males, and SL and SVL 
in females, respectively (Supplementary Table S2). The CAN2 
(explaining 12.1% and 22.9% of the whole variance in males and 
females, respectively) differentiated O. schmackeri and O. sp2, 
largely based on HW and FL in males, and NSD and FTL in 
females (Supplementary Table $2). 

The CDA, utilizing 18 morphometric variables from 
SMF6241(male) and the above male samples, was also 
conducted to reevaluate the taxonomic attribution of SMF6241. 


The discriminant functions yielded significant Wilks’ A (P < 
0.001 for two axes; Supplementary Table S3), and the accurate 
classification rate was 98.8% in raw data and 95.2% in leave- 
one-out cross-validation, respectively. The plot of SMF6241 was 
grouped with O. schmackeri and greatly varied from the other 
two groups (Figure 2C). The squared Mahalanobis distances 
from SMF6241 to О. schmackeri, to О. sp2, and to О. hejiangensis 
were 10.209, 30.809 and 97.481, respectively. 

Furthermore, we used univariate assessments to explore the 
morphological differences between O. sp2, O. schmackeri and O. 
hejiangensis (Table 2). The morphological descriptive statistics 
for each measurement are shown in Supplementary Table S4. 
For males, O. sp2 was extremely different from O. hejiangensis 
for all the 24 characters with nearly all P < 0.05 in the ANOVA, 
whereas only nine traits (SVL, HW, SL, NEL, ТЕЕ, TD, HLL, TL 
and TW) between O. sp2 and O. schmackeri were significantly 
different. For females, however, less morphological variance 
among species was present than in males, and it varied 
among species pairs. The O. sp2 can be distinguished from O. 
hejiangensis in 10 characters (SVL, HW, IAE, TED, LAHL, HAL, 
HLL, TL, FTL and FL) with P < 0.05, while it diverged from 
О. schmackeri in 8 traits ТОР, ТЕЕ, ТАЕ, HAL, HLL, FTL, FL 
and OMTL). Between sexes, the difference was significant in 
all the three species with females being larger than males in 
both body size and body shape (ANOVA, P < 0.001; Figure 1; 
Supplementary Table S5). However, O. sp2 exhibited distinct 
sexual dimorphism both in magnitude and direction. Although 
О. sp2 had the largest ratio (19) of female SVL to male SVL than 
О. schmackeri (18) and О. hejiangensis (1.4), it showed the fewest 
size-corrected traits which were significantly different between 
sexes than the other two species. Only seven traits (ie, IOD, 
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LAHL, HAL, HLL, TL, FTL and FL) in females were found 
to be significantly larger than in males in O. sp2 (ANOVA, 
P « 0.05; Supplementary Table S5), whereas more traits were 
significantly different between sexes in O. schmackeri (ie, HW, 
TD, TED, LAHL, HAL, HLL, TL, FTL, FL and OMTL) and O. 
hejiangensis (ie, HW, IOD, ED, TED, LAD, ГАНІ, HAL, HLL, 
TL, FTL, FL, TW and МТМ). In addition, О. sp2 displayed 
reversals of sexual dimorphism in head shape and lower arm 
diameter, but the differences were not statistically significant 
(Supplementary Table S5). 


32. Genetic structure with microsatellite DNA According 
to AK, the best-supported number of clusters in STRUCTURE 
was three (AK = 3; Figure 3). As shown in Figure 3, the entire 
population was divided into three genetic clusters. Nearly all the 
individuals of О, schmackeri (М = 48), O. hejiangensis (N = 21) and 
О. sp2 (М = 23) were assigned to one of the three homogenous 
gene pools (4 > 095; 1 = pure cluster). Only one individual of О. 
sp2 showed a slightly lower q value of 0.94. 


3.3. Vocalization Call parameter statistics are shown in Table 
3. We recorded and analyzed the spontaneous vocalizations 
from five males of O. sp2, two males of O. schmackeri and one 
male of O. hejiangensis, respectively. Male advertisement calls 
of these species were short and discontinuously separated 
by long silent intervals and were reminiscent of a bird's call 
"giu, qiu". The oscillograms and spectrograms of calls showed 
species-specific patterns (Figure 4). Male calls of O. sp2 and 
O. schmackeri were composed of 1-2 short single-pulsed notes 
with 4 frequency bands. In contrast, O. hejiangensis produced 
both short calls (Figure 4) and long calls (not shown here). The 
short calls of O. hejiangensis included 3-4 pulses with energy in 
2 frequency bands, and the long calls were composed of pulse 
series. 


In acoustic features, O. sp2 was distinct from its congeners 
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(Table 3). The average note duration of O. sp2 (mean 0.053 + 
0.018 s, rang 0.017-0.100 $, N = 77) was prominent shorter than 
both of О. schmackeri (mean 0.078 + 0.145 s, rang 0.053-0.115 N 
= 21) and O. hejiangensis (mean 0.121 + 0.016 s for short calls, rang 
0.102-0.145 s, N = 12; Kruskal-Wallis test, both P < 0.001). The 
second harmonic was the dominant frequency (DF) in all the 
three species. The DF of О. sp2 (3.679 + 0.460 kHz, М = 77) was 
significantly lower than O. schmackeri (3.965 + 0.272 kHz, N = 
21; Kruskal-Wallis test, P = 0.009) but much higher than that of 
O. hejiangensis (2.807 + 0.045 kHz, N = 12; Kruskal-Wallis test, P 
< 0.001). Significant differences in fundamental frequency (FF) 
were found only between О. sp2 and О. schmackeri (Р = 0.001). 
In addition, the average maximum frequency (MAXF) of О. 
sp2 (7.806 + 0.360 kHz, ranging 7.231-8913 kHz , М = 30 ) was 
lower than that of О. schmackeri (8.474 + 0.509 kHz, ranging 
7.786-9.336 kHz, М = 16; Kruskal-Wallis test, P = 0.003), but 
much higher than that of O. hejiangensis (Kruskal-Wallis test, P 
< 0.001). 

In brief, integrative evidence of morphology, molecules 
(mitochondrial and nuclear makers) and acoustics suggested the 
species status of O. sp2, and we proposed and described O. sp2 as 
Odorrana ichangensis Chen, sp. nov. 


3.4. Species Accounts 

Odorrana ichangensis Chen, sp. nov. (Figures 1A, 5-7; Table 4) 
Holotype HNNU18071139, an adult male (Figures 1A; 5A, В, 
E), collected by H. J. SHEN, M. Y. XU and C. KE on 24 July 
2018 from Gaojiayan Town (N30'3845.6', Е110°5777.2"; 323 m 
elevation), Changyang Co, Yichang City, Hubei Prov., China. 
Allotype HNNU18071155, an adult female (Figures 1A; 5C, D, 
Е), collected simultaneously and syntopically with the holotype 
(HNNU18071139). 

Paratypes Fifty male and thirty female adult specimens 
were collected by X. H. CHEN, L. QIAO, Y. Q. LU, Y. J. ZHU, 
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Table 2 Results of one-way ANOVA with corrected P values for differences in morphometric variables among O. schmackeri, O. hejiangensis 
and O. sp2, separated by sex. 


Males Females 
Variables 
SP vs. SM SP vs. HE SM vs. HE SP vs. SM SP vs. HE SM vs. HE 

SVL 0.0014** AON EO S 0.802 0.005** 0.002** 
HL 0.9067 <0.0001**** <0.0001**** 0.295 293 1 
HW =0000 ес <0.0001**** Eo 0.156 0.017* 0.401 
SL 0.0004*** <0.0001**** <0.0001**** 1 1 0.109 
NSD 0.411 <0.0001**** 0.0003*** 1 1 0.902 
NEL 0.005** <0.0001**** <0.0001**** 1 0.761 0.717 
IND 0.107 20.000122 <0.0001**** 0235 0.279 1 
IOD 0.432 <0.0001**** <0.0001**** 0.002** 0.331 0.008** 
IFE 0.015* <0.0001**** =0000 55 0.028* 0.254 0.403 
IAE 1 <0.0001**** <0.0001**** 0.049* 0.004** 0.367 
ED 0.541 < ООО <0.0001**** 1 0.239 0.055 
UEW 0.738 0.029* 0.0001*** 1 0.998 1 
TD 0.001** 0.0004*** 1 1 1 1 
TED 1 0.027* 0.004** 1 0.001** <0.0001**** 
LAHL 0.09 <0.0001**** 000012 0.083 0.004** 0.169 
LAD 0.278 <0.0001**** <0.0001**** 0.437 0.376 1 
HAL 0.745 <0.0001**** <0.0001**** 0.004** <0.0001**** 0.132 
HLL 0.023* <0.0001**** <0.0001**** 0.015* <0.0001**** 0.009** 
ШЕ, 00035 <0.0001**** =< (OOO Messrs 0.06 0.0002*** 0.005** 
TW <0.0001**** <0.0001**** <0.0001**** 0.228 1 0.028* 
FTL 1 <0.0001**** <0.0001**** 0.011* <0.0001**** 0.0003*** 
FL 1 <0.0001**** <0.0001**** 0.002** 0.001** 1 
OMTL 0.129 0.0001*** <0.0001**** 0.04* 0.182 0.907 
OMTW 0.165 <0.0001**** <0.0001**** 0.165 0.139 1 
Note: SM, O. schmackeri; HE, O. hejiangensis; SP, O. sp2. 
* P 0105; 
l 
n: р 0.001; 
****- P < 0.0001. 
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Figure 3 Bayesian clustering analysis based on ten microsatellite loci using STRUCTURE. Each individual is represented by а 
vertical bar and grouped by population and species (A). Different species are denoted as different colors (Purple: O. hejiangensis; 
Green: O. schmackeri; Red: O. sp2.). The optimal K value is three (B). 
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Q. Q. FENG, L. M. CHANG, Y. K. WAN, H. J. SHEN, M. 
Y. XU and C. KE from 2009 to 2018 from eight localities in 
China (Supplementary Table S6). АП the type specimens were 
deposited in the Herpetological Research Laboratory, College of 
Life Science, Henan Normal University. 

Diagnosis Odorrana ichangensis sp. nov. differs from its 
congeners by the following combination characters: (1) snout- 
vent length (SVL) of adult females approximately twice as long 
as males, 82.4 + 45 (72.4-930 mm, М = 30) in females, 429 + 24 
(38.2-471 mm, N = 50) in males; (2) head length greater than 
width in males, smaller than width in females; (3) tympanum 
large, greater than one-half diameter of the eye; (4) eyes big and 
prominent, width of upper eyelid (UEW) approximately 3/4 
of eye diameter (ED); (5) dorsolateral fold absent; (6) legs long, 
tibiotarsal articulation reaching to nostril when leg stretched 
forward; (7) dorsum emerald green to yellow green in life, 
scattered irregular dark brownish spots; (8) lips with black 
spots; (9) dorsum of limbs with black-brown transverse bars; 
(10) disks on all digits with lateroventral grooves; (11) relative 
finger lengths: III>IV>II>I; (12) feet fully webbed; (13) males 
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with paired subgular vocal sacs located at corners of throat; 
(14) pinkish lineae musculinae on dorsal side in males; (15) gray- 
white pads on dorsal surface of first finger in males; (16) white 
spinules on ventral surface of throat, chest, abdomen, limbs and 
upper and lower lips in males and females; (17) females having 
creamy yellow eggs without pigmentation, mean diameter of 
23 mm (Figure 5F). 
Description of holotype Body miniaturized (SVL 38.5 mm); 
head length equal to width; snout obtusely pointed in profile, 
projecting beyond lower jaw; nostril lateral, distance to snout tip 
equal to eye; canthus rostralis distinct; loreal region concave and 
oblique; internasal distance larger than interorbital distance; 
tympanum-eye distance small (TED/TD=0.3); tympanum 
distinct, rounded, 75% eye diameter, not depressed relative to 
skin of temporal region; tympanic rim elevated relative to 
tympanum; vomerine teeth well developed, on two oblique 
ridges, separated by a distance barely equal to length of one; 
tongue narrowly cordiform, deeply notched posteriorly. 

Hand moderately long, hand length (HAL) 03 times of SVL; 
tips of fingers moderately dilated, becoming pointed, cordiform 


8000 
B 
6000 
4000 
чес 
2000 
0 
0 0.1 0.2 0.3 
Time(s) 
8000 
D 
6000 4 
4000 4 
вв... 
2000 4 
0 т 
0 0.1 0.2 0.3 0.4 0.5 
Time(s) 
8000 
F 
6000 
t 
4000 ф 
а 
= 
2000 t 
0 
0 0.1 0.2 0.3 
Time(s) 


Figure 4 Oscillograms (left) and spectrograms (right) of the advertisement calls from О. schmackeri (A, B; air temperature = 26.2 °C; 
air humidity = 99.9%), O. hejiangensis (C, D; air temperature = 25.9 °C; air humidity = 94.2%) and O. sp2 (E, F; air temperature=26.2 


°C; air humidity = 96.4%). 


Table3 Acoustic parameters and P values in Kruskal-Wallis ANOVA on ranks adjusted by Bonferroni-Holm for repeated comparison of call 
traits among О. schmackeri, О. hejiangensis and О. sp2. Values are presented as mean + standard deviation (SD) followed by range. 


P values 
SM vs. SP HE vs. SP SM vs. HE 


0.078 + 0.145 0.121 + 0.016 0.053 + 0.018 
Note duration (s) (0.053-0.111) (0.102-0.145) (0.017-0.100) <0.001 <0.001 0.067 
(N = 21) (N = 12) (М 77) 


2.152 + 0.150 1.800 + 0.030 1.747 + 0.355 
FF (kHz) (1.985-2.411) (1.756-1.815) (1.216-2.617) 0.001 1 0.335 
(N = 8) (№ = 4) (N = 40) 


Note: Bold values indicate significant effects. SM, HE and SP represent O. schmackeri, O. hejiangensis and O. sp2, respectively. DF and FF 
indicate dominant frequency and fundamental frequency, respectively. MAXF represents the maximum frequency. 


Table 4 Morphometric measurements (in mm) of the type series of O. ichangensis sp. nov. from eight populations across its distribution: mean 
+ SD (range) and ratio of each measurement to the SVL (%). 
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disks with lateroventral grooves on all fingers; width of finger 
III disc slightly larger than width of distal phalanx (FPW/ 
ОРМ/=1.1), approximately 55% the dia meter of tympanum; 
fingers bearing narrow lateral fringes; relative finger length: 
ШЪ>ТУ>П>!, subarticular tubercles ovoid, markedly elevated and 
prominent: 1, 1, 2, 2; weak supernumerary tubercles below base 
of finger II, III and IV; inner metacarpal tubercle oval, elongated, 
larger than length of finger I; outer two metacarpal tubercle 
small, feebly prominent; fingers without web (Figures 5A, B). 

Hind limbs comparatively robust and long, tibia length 63% 
of SVL; tibiotarsal articulation of the adpressed limb reaching to 
nostril; foot slightly shorter than tibia (FL/TL-0.9); relative toe 
lengths: IV» V2IIDIDT, tip of toe V extending to the distal edge 
of distal subarticular tubercle on toe IV; tarsal fold absent; tips 
of toes tapering dilated into narrow and long triangular disks; 
lateroventral grooves present; feet entirely webbed to disks; web 
becoming narrower and reaching the base of disk except toe IV 
webbed to distal subarticular tubercle; subarticular tubercles 
prominent, ovoid; inner metatarsal tubercle elongated; outer 
metatarsal tubercle absent; inner tarsal fold absent (Figures 
5A, B). 

Skin on dorsum slightly shagreened; dense granules on the 
entire dorsum and upper sides of limbs; small round tubercles 
on lores, temporal region, and posterior to tympanum; larger 
flat tubercles on flanks of trunk from back of eye to base of 
forelimb; no dorsolateral fold; two rictal glands; small white 
spinules on ventral surface of throat, chest, abdomen, limbs and 
upper and lower lips. 

Color of holotype in life Dorsum emerald green; small 
and irregular dark brownish spots on dorsal surface of head, 
back, limbs and flank; small and black spots on the margin of 
snout tip and lower jaw; dorsal surface of limbs emerald green 
with narrow, diffuse dark brown crossbars (4—5 strips on 
thigh), interspersed with small dark brown spots; tympanum 
brownish, center and margin yellow; pinkish lineae musculinae 
on dorsal side; venter ivory to yellowish-brown with diffuse 
brown markings on chest; brownish pigmentation on gular 
pouch. 

Color of holotype in preservative Dorsum black-brown; side 
of head, supratympanic fold, upper and lower lip yellowish- 
gray with small dark brown spots; flanks whitish-gray with 
fewer light brown spots and creamy-white flat tubercles; rictal 
glands creamy-white; nuptial pad whitish-gray; tympanum 
brown, margin dark brown; dorsal limbs with black-brown 
spots and banding; venter yellowish-gray with diffuse dark 
brown markings on gular pouch and chest. 

Sexual dimorphism Males small (SVL, mean 42.9 + 24 mm, 
range 38.2—47.1 mm), on average half of females (SVL, mean 
82.4 + 45 mm, range 724-930 mm), but obviously larger than 
females in the relative size of ED and TD (Table 4); nuptial pads 
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whitish-gray develop on the inner aspect of the first fingers, 
paired subgular vocal sacs located at corners of throat, pinkish 
lineae musculinae on dorsal side in males. 

Variation in the type series Variation in morphometric 
measurements of the type series from eight populations across 
its distribution range are listed in Table 4. Basically, the type 
series agrees with the holotype coloration, but also shows 
differences in color pattern and spots (Figures 6-7). The venter 
surface is grey white with light brown mottling on the lower 
jaw, throat, chest and upper abdomen in male HNNUI0071063 
(Figure 6F). The color of the dorsum surface is gray-brown in 
female HNNU1408II140 (Figure 7G). Large dark brown spots 
are present on the dorsum of the posterior body in female 
HNNUI5071304 (Figure 7D), and the lateral side of flanks 
formed networks in female HNNUI4081II133 (Figure 7F). 
Etymology The specific epithet “ichangensis” refers to the type 
locality of the new species as well as O. schmackeri, Yichang City, 
Hubei Province, China. The "ichang" is an old spelling and a 
transliteration for “Yichang”. The suggested common English 
name is the "Yichang Odorous Frog" for this species. 
Distribution and ecology Odorrana ichangensis sp. nov. 
inhabits in mountain streams in western China, distributed in 
Qionglai Mountain, north of Dalou Mountain, Wushan and 
Daba Mountains, from the western of Sichuan Basin along its 
southern edge up to the northeastern edge, reaching the eastern 
edge of Wuling Mountain. 

АП specimens were collected at night on rocks and 
stream banks along slow flowing canyon streams in broad- 
leaved forests at 246-1025 m elevations (Figure 8). Breeding 
season was approximately from late July to Mid-August. The 
amplexus and spawn behaviors were found on July 27, 2018 at 
Changyang Со, Hubei Prov. (N30°3623.97’, E1112:588', 260-336 
m elevations; Figure 1A), July 29, 2015 at Xuyong Co, Sichuan 
Prov. (662-747 m elevations), and August 16, 2014 at Emei 
Mountain, Sichuan Province of China (611-672 m elevations). 
Adult female spawn 1300-1784 eggs once a time with a mean 
diameter of 229 mm. Mature eggs are creamy yellow, without 
pigment (Figure 5F). 

Comparisons Odorrana ichangensis sp. nov. closely resembles 
O. schmackeri in morphology, and both were collected from 
Gaojiayan Town, Hubei Province. simultaneously and 
syntopically. However, O. ichangensis sp. nov. differs from O. 
schmackeri by having the dorsum green with small, irregular, 
dark brown spots (dorsum green with distinct, round, large 
light brown spots in O. schmackeri), having males with greater 
relative head length, head width and internarial distance than 
females, having supernumerary tubercles and lateral fringes 
on toe V (absent in O. schmackeri), absence of stripes on canthus 
rostralis (O. schmackeri with narrow brown stripes on canthus 
rostralis, from tip of snout along margin of upper eyelid to 
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Figure 5 (A, B) Dorsolateral and ventral views of the holotype 
(HNNU18071139, male). (C, D) Dorsolateral and ventral views 
of the allotype (HNNU18071155, female). (E) Dorsolateral 
view of the amplexus of the holotype (upper) and allotype 
(lower). (F) Eggs of female O. ichangensis sp. nov. Photographs 
by H. J. SHEN. 


posterior corner of upper eyelid), the absence of skin fold 
(present in O. schmackeri), and females with pigmentless eggs 
(eggs beige with black poles in O. schmackeri). 

Odorrana ichangensis sp. nov. differs from O. hejiangensis by 
having a greater ratio of female SVL to male SVL (19 in O. 
ichangensis sp. nov., 14 in О. hejiangensis, males SVL 473-57.8 
mm, mean 52.7 mm, females 679-815 mm, mean 759 mm), 
having males with greater relative head length, head width 
and internarial distance than females, having the dorsum green 
with narrow crossbars on limbs (anterior half of dorsum green, 
posterior half of dorsum brown with wide crossbars on limbs 
in O. hejiangensis), having small, irregular, dark brown spots on 
dorsum (distinct, round, large brown spots on dorsum, blending 
into shadow at rear of back), the absence of white spinules on 
dorsum (О. hejiangensis with white spinules on surface of upper 
eyelid, limbs, tympanic region and rear of back; Deng and Yu, 
1992). 

Odorrana ichangensis sp. nov. differs from O. kweichowensis 
(Li et al., 2018) by the longer snout-vent length of females 
(O. ichangensis sp. nov. 72.4-93.0 mm, mean 82.4 mm; O. 
kweichowensis females SVL 62.4-81.1 mm, mean 73.6 mm), 
having three outer metacarpal tubercles (two outer metacarpal 
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Figure 6 Variation of males in morphology and color 
patterns of O. ichangensis sp. nov. in life. (A, B) Dorsolateral 
views of HNNU15071067 and HNNU15071061 from Gufu 


Town, Xingshan County (Co.), Hubei Province (Prov.). (C) 
Dorsolateral view of HNNU15071104 from Mawu Town, 
Shizhu County (Co.), Chongqing City. (D) Dorsolateral view 
of HNNU15071275 from Shuiwei Town, Xuyong County 
(Co.), Sichuan Province (Prov.). (E, F) Dorsolateral view and 
ventral view of HNNU10071063 respectively from Gaojiayan 
Town, Changyang County (Co.), Hubei Province (Prov.). 
(С) Dorsolateral view of HNNU1608II239 from Wenquan 
Town, Kaizhou District, Chongqing City. (H) Ventral view of 
HNNU1408I1162 from Bifengxia Town, Yaan City, Sichuan 
Province (Prov.). Photographs by X. H. CHEN, L. QIAO, Y. J. 
ZHU and Q. Q. FENG. 


tubercles in O. kweichowensis), the second finger longer than the 
first (the second finger shorter than the first in O. kweichowensis), 
lower arm and hand shorter than a half SVL (lower arm and 
hand longer than a half of SVL in O. kweichowensis), tibiotarsal 
articulation of adpressed limb reaching nostril level (tibiotarsal 
articulation of adpressed limb reaching to the level between eye 
to nostril in О. kweichowensis). 

Odorrana ichangensis sp. nov. differs from O. huanggangensis 
(Chen et al, 20102) by having the dorsum green with small, 
irregular, dark brown spots (dorsum yellow-green with distinct, 
round, large brown spots in O. huanggangensis), toe IV webbed 
to distal subarticular tubercle (all toes webbed to base of disc 
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Figure 7 Variation of females in morphology and color 
patterns of O. ichangensis sp. nov. in life. (A) Dorsolateral view 
of HNNUI16081I147 from Jiucheng Town, Daozhen County 
(Co.), Guizhou Province (Prov.). (B) Dorsolateral view of 
HNNU140711126 from Wenquan Town, Kai County (Co.), 
Chongging City. (C) Dorsolateral view of HNNU15071063 
from Gufu Town, Xingshan County (Co.), Hubei Province 
(Prov.). (D, E) Dorsolateral views of HNNU15071304 and 
HNNU15071276 respectively from Shuiwei Town, Xuyong 
County (Co.), Sichuan Province (Prov.) (F-H) Dorsolateral 
views of HNNUI14081I133 and HNNU1408II140 and ventral 
view of HNNUI140811155 respectively from Bifengxia Town, 
Yaan City, Sichuan Province (Prov.), China. Photographs by X. 
Н. CHEN, Г. QIAO, Y. J. ZHU and О. О. FENG. 


in O. huanggangensis), females with a larger ratio of ED/SL (O. 
ichangensis sp. nov., 0.58-0.86, mean 074; O. huanggangensis, 
0.52-0.63, mean 0.58), females with a lower ratio of TED/SVL 
(О. ichangensis sp. nov., 0.03-0.05, mean 004; О. huanggangensis, 
mean 0.05), females with pigmentless eggs (eggs beige with 
black poles in O. huanggangensis) and the absence of stripes on 
canthus rostralis (narrow brown stripes on canthus rostralis, 
from tip of snout along margin of upper eyelid to posterior 
corner of upper eyelid in O. huanggangensis). 

Odorrana ichangensis sp. nov. differs from O. tianmuii (Chen 
et al, 2010b) by having the dorsum green with small, irregular, 
dark brown spots (dorsum yellow-green with distinct, brown 


spots in О. Наптий), tibiotarsal articulation of adpressed limb 


Figure 8 Habitats of O. ichangensis sp. nov. (A) Habitats 
in the type locality, Gaojiayan Town, Changyang County 
(Co.), Hubei Province (Prov.). (B) Habitats in Gufu Town, 
Xingshan County (Co.), Hubei Province (Prov.). (C) Habitats 
in Wenquan Town, Kai County (Co.), Chongqing City. 
(D) Habitats in Mawu Town, Shizhu County (Co.), 
Chongqing City. (E)Habitats in Jiucheng Town, Daozhen 
County (Co.), Guizhou Province (Prov.). (F) Habitats in 
Shuiwei Town, Xuyong County (Co.), Sichuan Province 
(Prov.). (G) Habitats in Huangwan Town, Emeishan City, 
Sichuan Province (Prov.). (H) Habitats in Bifengxia Town, 
Yaan City (Со.), Sichuan Province (Prov.), China. Photographs 
by X. H. CHEN, L. QIAO, Y. J. ZHU and Q. Q. FENG. 


reaching nostril level (tibiotarsal articulation of adpressed 
limb reaching to the level between eye to nostril in O. tianmuii), 
females with lower ratio of TD/ED (О. ichangensis sp. nov., 
0.42-0.64, mean 0.52; О. tianmuii, 0,58-0.78, mean 0.63), females 
with pigmentless eggs (eggs beige with black poles in O. tianmuii) 
and the absence of stripes on canthus rostralis (narrow brown 
stripes on canthus rostralis, from tip of snout along margin of 
upper eyelid to posterior corner of upper eyelid in O. tianmuii). 


4. Discussion 
4.1. Taxonomic confirmation and species delimitation The 


genus Odorrana has a very complicated history of systematics, 
and this is partly due to insufficient type series or unreliable 
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published host records. For example, Rana nebulosa was erected 
from a juvenile specimen and the holotype specimen was 
lost (Hallowell, 1861; Bain et al., 2003), O. sinica lacked definite 
type locality and detailed description (АМ, 1927 “1925”) and 
O. livida used to be short of description of male specimens 
until Rujirawan et al. (2018) published a detailed description 
(Rujirawan et al., 2018, Chen et al., 2020). Odorrana ichangensis sp. 
nov. most closely resembles O. schmackeri, and they co-distribute 
with O. hejiangensis at the type locality of O. schmackeri 
(Gaojiayan Town, Hubei Province). Therefore, it is essential 
to reevaluate the species belonging of the SMF6241 before we 
describe the new species from the same type locality. Here, we 
uncovered the fine-scale divergence between O. ichangensis sp. 
nov. and its sympatric species O. schmackeri and O. hejiangensis 
from Gaojiayan with multiple data sets. 

Our results revealed marked differentiation in morphology 
between O. ichangensis sp. nov. and its sympatric congeners. 
For example, O. ichangensis sp. nov. was clearly separated from 
the other two species in the CDA. In addition, O. ichangensis 
sp. nov. showed the opposite direction in sexual dimorphism 
in head shape and forelimb size, when compared to the other 
two species. Males of O. ichangensis sp. nov. seemed to have a 
relatively (size-corrected) larger head (e.g, HL, HW and IND; 
Table S4, S5) than females, and the size-corrected LAD (ie, 
lower arm diameter) of females was unexpectedly larger than 
that of males. Bayesian clustering method (STRUCTURE) 
using microsatellite data further emphasized three different 
genetic clusters. In addition, O. ichangensis sp. nov. also showed 
species-specific characteristics based on the bioacoustic analyses. 
As call signals are considered to be an important premating 
isolation mechanism for frogs (Leary, 2001; Wells, 2007), 
and STRUCTURE detected no heterospecific gene flow, 
we can expect that there is a complete reproductive barrier 
between ichangensis sp. nov. and its close relatives. Thus, the 
divergences revealed based on the nuclear gene, bioacoustics 
and morphology data support the validity of the new species 
O. ichangensis sp. nov. This work represents an example 
to resolve the taxonomy of old museum specimens under 
historical taxonomic conflicts. It also highlights the need for the 
integrative framework to address species boundaries, especially 
when morphological evolution is conservative. 


42. The association between morphological and ecological 
features Our results support the idea that ecological similarity 
weighed more with trait convergence than phylogenetic 
relatedness (Liu et al., 2013; Jiang et al., 2018). Previous 
phylogenetic studies have revealed that O. ichangensis sp. nov. 
was more closely related to O. hejiangensis (with a Kimura 
2-parameter distance of 3.996) than to O. schmackeri (with a 
Kimura 2-parameter distance of 6.8%; Chen et al., 2013; Zhu, 
2016). However, O. ichangensis sp. nov. closely resembles O. 
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schmackeri than O. hejiangensis in morphometrics and color 
patterns (Figure 1). The convergent evolution of phenotypes 
has been identified associated with similar environments in 
amphibians (Moen et al., 2016) and reptiles (eg. Anolis lizards, 
Losos et al., 1998 and Mahler et al., 2013, and Cryptoblepharus 
lizards, Blom et al., 2016). Odorrana hejiangensis always inhabits 
on mountain rocky substrates such as cliffs and huge rocks near 
rushes and narrow rivers in closed and moist microhabitats 
(Deng and Yu, 1992), while O. ichangensis sp. nov. and O. 
schmackeri (Li et al, 2017), showing similar habitat preferences, 
dwell on barer stones or rocks near slow and wide rivers in 
open and more arid environments (Figures 1, 8). The pattern 
of interspecific variance in traits may suggest the important 
roles of divergent selective pressures in driving trait evolution 
in this group. Previous studies (Mcguigan et al., 2005; Conte et 
al., 2012; Moen et al., 2016; Laiolo, 2017) suggested that similar 
environmental habitats could lead to phenotypical similarity 
among independently evolutionary species due to the same 
selective pressures, while heterogeneous ecological conditions 
might promote genetic and phenotypic diversification among 
lineages. It is likely that identical selective pressure caused 
by the same environmental factors promoted morphological 
convergence in O. ichangensis sp. nov. and O. schmackeri. In 
contrast, divergent selection pressures from heterogeneous 
environments might drive morphological divergence in O. 
hejiangensis. Such environment-driven interspecific phenotypic 
differentiation has been revealed in Australasian warblers 
(Acanthizidae; Garcia-Navas et al., 2018) and African drongos 
(Dicrurus adsimilis; Fuchs et al., 2018). 


4.3. Ecological adaptation of signals Recent studies have 
suggested that signal ecology might play an important role 
in the diversification of signal traits (Kolm et al., 2012; Snell- 
Rood, 2012; Odendaal et al., 2014; White and Kemp, 2015; 
Maan et al., 2017; Maluleke et al., 2017), through a process 
termed sensory drive (Endler, 1992; Endler and Basolo, 1998). 
Under the scenario of sensory drive, environmental factors 
(e.g. temperature and humidity) can exert selective pressure 
on signals as mediated through sound absorption (Fuller 
and Behavior, 2010; Malone et al., 2014). In particular, signals 
of lower frequency and longer duration could benefit the 
animals in humid environments (Snell-Rood, 2012; Jiang е al., 
2015; Maluleke et al., 2017). Here, the revealed species-specific 
signal features agreed with this hypothesis. In the field survey, 
O. ichangensis sp. nov. and O. schmackeri were always found 
in open (meaning higher temperature) and relatively arid 
microhabitats, while О. hejiangensis occurred in closed (meaning 
lower temperature) and moister habitats. Correspondingly, O. 
ichangensis sp. nov. and O. schmackeri emitted calls with higher 
frequency (the average DFs were 3.679 kHz and 3.965 kHz, 
respectively) and shorter duration (the average note duration 
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was 0.053 s and 0.078 s , respectively), whereas O. hejiangensis 
produced calls with lower frequency (DF 2.807 kHz) and longer 
note duration (0.121 s). In addition, O. hejiangensis showed a more 
complicated call structure. For example, calls of O. hejiangensis 
were 3- or 4-pulsed and showed both short calls and long calls 
(approximately 5 s per call). 

On the other hand, the divergence in temporal and spectral 
properties of male calls might result from ecologically associated 
traits, such as body size (Ophir et al., 2010; Snell-Rood, 2012; 
Odendaal et al., 2014; Lee et al., 2016; Maluleke et al., 2017). 
Interestingly, we found significant interspecific differentiation 
of body size in males but no differences in females. Male O. 
hejiangensis (mean SVL 52.7 mm) is obviously larger than male 
O. ichangensis sp. nov. and O. schmackeri (mean SVL is 42.4 
mm and 446 mm respectively). As a proxy of environmental 
quality, body size in relation to many other ecological traits 
always bears strong selection and, by extension, might drive 
the divergence in signals by imposing effluence on signal traits 
(eg, dominant frequency and call duration, Stewart et al, 2016; 
Köhler et al., 2017; Derryberry et al., 2018). In general, the call 
duration correlated positively with body size (Castellano and 
Cristina, 1998), whereas the dominant frequency correlated 
negatively (Hoskin et al., 2005; Pfennig and Pfennig, 2005; Nali 
and Prado, 2014). In other words, the heterogeneousness of the 
microenvironment could also indirectly lead to interspecific 
variation in acoustic features via association with body size. 
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Appendix 


54. Вапа schmackeri Bttgr. n. sp.” 


1054,2a. Erw. Kao-cha-hien hei Iehang, Prov. Hubei, Contr-China. Gesch, 
1991 v. В. Schmackor, Shanghai 


35. Rana leptodaetyla Bigr. SME (ro 
Bonlenger, Cat. pag. 57, ~ 
1057 а. Зе, Kodaikanal, Palni Hills, S.-Imlien. (esch. 1889 v. Theodor 
Kolb. Madras, 


е) Char. Plica lateralis glandalosa nulla. apices digitornm manns et 
pedis distincte. dilatati, diseis parvis. Lingua поп papillifera. Турат 
distinetissimum. / magnitutidinis bulbi. Cutis natatoria fere completa, Series 
dentium vomeris perobliqui, valde inter se separati. medias inter ehoanas hand 
procul ab angulo anteriore et interiore choanarum incipientes, Cutis laevis. 

Vomerziihne in zwei kräftigen, sehr schiefen Reihen, die zwischen den 
inneren Nasenüfinungen nicht weit von ihrem oberen Innenrande beginnen 
nnd erheblich hinter das Nivean derselhen nach hinten reichen. Kopf von 
mäßiger Grüße, doch im der Ohrgegend fast doppelt ко breit wie der Rumpf 
in der Krenzbeingegend ; Schnauze kanm länger als die Orbita, leicht zugespitzt 
mit stumpfer Zügelkante und schief al fallen lor. etwas ausgehöhlter Ziigelgexend ; 
Nasloch der Schnanzenspitze etwas näher ts dem Auge: Jnterorbitalrawm wenig 
schmäler als ein einzelnes Augenlid. Trommelfell sehr dentlich, von „ Angen- 
größe. Finger mäßig schlank, erster etwas linger als der zweite; Zehen 
mäßig lang, Schwimmhant vollkommen, nur das letzte Glied der vierten Zche 
herausragend, aber beiderseits gesimnt; Finger- und Zehenspitzen in kleine — 
Haftscheiben verbreitert; Subartienlarknötchen mäßig groß; ein einziger 
innerer, schr schwacher nnd undeutlicher Metatarsaltuberkel. Tarsalfalte schwach 
angedentet. Das Hinterbein reicht, nach vorn gelegt, mit dem Tihiotarsakrelenk 
bis zur Schnanzenspitze; die Tilia besitzt erheblich mehr als die halbe Кор!- 
rumpflänge. Haut glatt oder fein lederartig gonarbt; eine schwache Drüsen- 
falte über dem Trommelfell, eine stärkere, weilgefürhte als Fortsetzung der 
Oberlippe unter dem Trommellell; hinterer Teil des Bauches mit schwacher, 
drüsiger Granulation. 

Kopf und Rumpf oben bläulich griünzran mit schwarzen Marmorzeich- 
nungen und grüßeren gerundeten Inselfierken, die zumeist etwas helier um- 
sänmt sind; Seiten oben gran, unten gelbweiß mit scharfbezrünzten schwarzen 
Rundflecken; Gliedmaßen braungran nnd nach den Seiten hin sehmntzizzclb 
mit sehr zahlreichen schwarzen, zumeist aus Flecken nnd Макеш bestehenden 
Querbinden übers&et, Kopfseiten hell graulichgriin mit großen schwarzen 
Flecken, die Drüsenfalte unter dem Trommelfell etwas lebhatter, zrünlichweiß 
gefärbt. Unterseite einfarbig schmutzigweiß, der Kinnraml init, undentlichen 
graulichen Wiirfelflecken ; die der Gliedmaßen schmutzig orangegelh. 

Maße. Kopfrumpflänge 42 mm  Ninterglielmaßen. . 75½ mm 


Kopfbreite . . . М, Oberschenkel . 22½ „ 
Augendnrehimesser б'/а „ Unterschenkel . . 96 b 
Trommelfell 4 , Breite der Scheibe des 

Vorderglielmaßen 28'/s „ 4. Fingers 1 p 


Figure 81 The original description of the holotype of O. schmackeri. 
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Male Female 
Variables 
CAN 1 CAN2 CAN 1 CAN2 
r (P) 266.031(<0.001) 75.915(<0.001) 114.103(<0.001) 40.359(0.014) 
Canonical correlation 0.967 0.815 0.939 0.83 
SVL 0.245 0.432 0.954 -0.095 


E 


-0.081 -1.29 -0.448 0.368 


2 
Ра 
о 


-0.172 0.051 0.542 -1.502 


— 


ND 0.039 -0.001 0.092 0.616 


— 
29| 
tri 


-0.396 -0.318 -0.518 -0.168 


т 
о 


-0.017 0.168 


о 
л 
© 


-0.677 


= 
J 


-0.636 -0.42 0.59 -0.332 


Е 


0.369 0.009 0.456 -0.604 


= 


0.104 0.159 -0.282 


© 
N 
N 


= 
ЕЗ 


0.699 -0.209 0.491 -0.32 


№] 


TL 0.156 -0.062 0.392 0.929 


OMTL 0.19 0.388 -0.144 0.655 


Note: Bold values indicate the largest two loading on each axe within each sex. 


Table $3 Standardized canonical coefficients from 24 morphometric characters, the Wilks Lambda, the canonical correlation, Ше eigenvalue 
and the percentage total variance. 


Variables CANI CAN2 


X (P) 258.273(«0.001) 67.712 («0.001) 


Canonical correlation 0.963 0.779 


I 
е 


-0.061 0.384 


n 
ES 


0.499 -0.416 


Z 


EL 0.085 0.044 


IOD -0.01 0.205 
unn ce „„ 
IAE -0.592 -0.279 
юв 0028 
UEW 0.298 0.086 
WW o:: 3 а 
ТЕР 0.624 -0.595 


© 


0.212 -0.072 


di 


LL 0.247 0.391 


-0.061 -1.105 


Е 


esl 
А 


0.499 0.112 


= 


MTW 0.085 0.016 


Note: Bold values indicate the largest two loading on each axe within each sex. 


Table S4 Measurements of 24 morphological characteristics of adult О. schmackeri, О. sp2 апа О. hejiangensis. 


ТИРУ О. schmackeri О. hejiangensis О. sp2 
$ (№31) 2 (N-23) $ (N=33) $ (N-21) $ (N-20) 2 (N-5) 
з 4464 801836 527426 759841 QAS20 824434 
40.4-47.9t 73.7-87.0 47.3-57.8 67.9-81.5 38.4-45.6 78.9-87.0 
н. їз EEE 2714818 0167409 fr 285426 
13.9-19.0 24.4-29.6 17.4-21.4 24.8-29.9 14.9-17.9 24.9-31.0 
EEUU оз» 37208 I 
HW 16.5 + 0.7 28.0 + 1.3 18.9 + 1.1 27.4 + 1.2 15.1 + 1.0 29.3 + 1.4 
маз 259-3L0 168-000 254-293 133-169 2738-3072 
37.00% 35.00% 36.00% 36.10% 35.60% 35.60% 
SL 72404 121404 84405 117405 6704 119407 
6.3-8.1 11.2-13.1 7.5-9.3 10.8-12.6 6.1-7.4 11.0-12.6 
ж lll l1 211 1 
NSD 3.24 0.4 5.14 0.5 3.7 + 0.5 5.0 + 0.4 3.0 +0.4 5.0 + 0.6 


7.20% 6.40% 7.00% 6.50% 7.20% 6.00% 


2.7-4.6 5.6-7.1 3.7-5.5 5.1-6.9 2.9-4.1 5.4-7.0 


2 
=, 


4.6 + 0.4 7.5 + 1.0 5.3 + 0.4 7.5 + 0.6 4.3 + 0.5 8.2 + 0.5 


10.30% 9.40% 10.10% 10.00% 10.20% 10.00% 


3.4-4.7 4.9-7.8 4.0-6.0 5.5-8.2 3.0-4.4 6.6-7.9 


= 


FE 7.8 + 0.7 12.9 + 1.1 9.1+0.7 13.4+0.9 7.3 + 0.7 14.3 + 1.3 


17.60% 16.20% 17.30% 17.60% 17.20% 17.40% 
IAE 1 %%% 140 —— 136210 18941 118209 2094140 
10.3-13.2 17.3-21.4 11.2-15.8 16.3--213 10.2-13.7 18.9-22.3 
ПО  26208 мш I 24808 EE M 
ЕР 6.4 + 0.5 9.8 + 0.8 7.3+0.7 9.2 + 0.7 6.2 + 0.6 9.9+1.0 
EM UU S вва 
14.4096 12.2096 13.8096 12.10% 14.60% 12.00% 
UEW 4404 68:05 EOS 6704 46205 70а | 
3.6-5.2 5.2-7.8 3.8-6.3 6.0-7.8 3.7-5.5 5.4-8.0 
EN 1 1 2 5 BOK 
TD 3.5 + 0.3 4.4 + 0.6 3.6 + 0.4 4.4 + 0.5 3.2 + 0.3 4.5 + 0.8 
м BAST 2842 3558 2734 3886 0 
7.90% 5.60% 6.70% 5.80% 7.50% 5.50% 
ТЕР 1.2 + 0.2 3.14 0.5 1.5 + 0.3 2.3 + 0.4 1.2 + 0.3 3.3 + 0.9 


2.70% 3.90% 2.80% 3.00% 2.90% 4.00% 


17.7-22.4 33.2-39.0 22.5-27.2 31.0-38.6 17.8-21.4 34.3-40.8 


Continued Table S4 


O. schmackeri O. hejiangensis O. sp2 


Traits „...... AAA ини ы Te 
$ (N-31) 2 (N-23) $ (N=33) $ (N-21) $ (М=20) 2 (N-5) 


3.0-5.0 4.4-7.7 4.0-6.1 5.2-6.8 3.1-4.6 6.0-7.01 


E 


14.0 + 0.7 23.2 + 0.9 16.3 + 0.5 22.6 + 1.0 13.8 + 0.6 25.0 + 1.4 


31.50% 29.00% 31.10% 29.80% 32.60% 30.40% 


72.2-86.8 134.2-154.7 89.9-105.1 131.3-153.5 73.45-84.69 148.4-158.1 


4 
ча 


25.5 + 1.1 47.0 + 1.6 30.7 + 0.9 45.3 + 1.5 24.6 + 0.9 49.2 + 3.1 


57.40% 58.70% 58.40% 59.80% 58.10% 59.80% 


4.8-6.7 8.9-11.9 5.6-7.5 9.4-12.2 3.4-5.7 10.0-12.2 


= 


33.9: 1.7 60.5 + 2.1 41.0 + 1.5 57.8 + 1.8 33.7 £1.5 63.7 + 3.0 


76.20% 75.60% 78.00% 76.30% 79.50% 77.60% 


21.1-26.3 41.0-46.3 26.7-31.8 39.8-47.0 22.9-26.7 43.4-50.2 


OMTL 2.6 + 0.3 4.3 + 0.4 3.2 + 0.3 4.4 + 0.5 2.8 + 0.3 4.9 + 0.5 


5.80% 5.40% 6.10% 5.90% 6.50% 5.90% 


0.7-1.2 1.5-2.2 0.9-1.6 1.5-2.3 0.6-1.3 1.7-2.7 


Note: t, ł and $ represent means + SD, range and ratio of each measurement to SVL respectively. Numbers in parenthesis next to symbols 6 
(males) or 9. (females) represent the number of specimens analyzed. 


Table S5 P values in Dunnett's post hoc test Юг one-way ANOVA and ANCOVA using SVL as covariates from О. schmackeri, O. sp2 and О. 
hejiangensis in separate sex. 


O. schmackeri O. hejiangensis O. sp2 
Raw data /SVL Raw data /SVL Raw data /SVL 


Traits 


HL <0.001Е 0.507 Е <0.001Е 0.086 Е <0.001Е 0.562 М 


n 
Е 


<0.001Е 0.339 Е <0.001Е 0.105 Е <0.001Е 0.269 Е 


2 


ЕТ, <0.001Е 0.758 М <0.001Е 0.112 Е <0.001Е 0.475 Е 


5 
=) 


<0.001Е 0.051 F <0.001Е 0.001 Е <0.001Е 0.045 Е 


> 
& 


«0.001F 0.580 F «0.001F 0.837M «0.001F 0.831 M 


= 


EW <0.001F 0.638 F <0.001F 0.531 М <0.001Е 0.845 М 


= 
т] 
Iv] 


«0.001F 0.004 F «0.001F 0.010 F «0.001F 0.879 Е 


= 
5 


<0.001F 0.056 M <0.001F 0.002 M <0.001F 0.204 F 


x 


LL <0.001F <0.001F <0.001F <0.001F <0.001F 0.001 F 


Е 


<0.001Е 0.654 М <0.001Е 0.065 Е <0.001Е 0.381 М 


= 
с 


<0.001Е <0.001Е <0.001Е <0.001Е <0.001Е 0.003 Е 


-- 


MTW «0.001F 0.311 F «0.001F 0.080 Е «0.001F 0.649 F 


Note: Bold values indicate significant effects (significant level is 0.05). Letters followed by values present the larger sex. F: females; M: males. 
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